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Abstract. The local thermodynamics of a system with long-range interactions 
in d dimensions is studied using the mean field approximation. Long-range 
interactions are introduced through pair interaction potentials that decay like a 
power law in the interparticle distance. We compute the local entropy, Helmholtz 
free energy and grand potential per particle in the microcanonical, canonical and 
grand canonical ensembles, respectively. The one-particle distribution function 
approach is also considered. From the local entropy per particle we obtain the 
local equation of state of the system by using the condition of local thermodynamic 
equilibrium. This equation of state turns out to be the ideal gas equation 
of state regardless of the pair interaction potential characterizing the long- 
range interactions. By volume integration of the relation between the different 
thermodynamic potentials at the local level, we find the corresponding equation 
satisfied by the potentials at the global level. It is shown that the potential energy 
enters as a thermodynamic variable which modifies the global thermodynamic 
potentials. As a result, we find a generalized Gibbs-Duhem equation which 
relates the potential energy to the intensive variables of the system. For the 
marginal case where the power of the decaying interaction potential is equal to 
the dimension of the space, the usual Gibbs-Duhem equation is recovered. As 
examples of application of this equation, we consider spatially uniform interaction 
potentials and the self-gravitating gas. We also point out a close relationship with 
the thermodynamics of small systems. 
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1. Introduction 

A great variety of systems in nature are dominated by long-range interactions. 
Examples are self-gravitating systems [iHl], two-dimensional vortices [5], nuclear 
physics [6J and also toy models such as the Hamiltonian mean field model [7]. Among 
other physical properties, the remarkable thermodynamic behaviour of systems with 
long-range interactions makes them extremely attractive. They are intrinsically non- 
additive and may have negative heat capacity in the microcanonical ensemble leading 
to ensemble inequivalence [5kl0|. For very instructive reviews on the subject with 
more examples and phenomena we refer to pT}llll| . 

Interactions in this kind of systems are characterized by slowly-decaying pair 
interaction potentials that couple the constituent parts of the system at large distances. 
Formally, a potential that decays a.sl/r'^ is said to be long-range iiv < d, where d is the 
dimension of the embedding space |10j . Systems having such interaction potentials are 
sometimes called systems with strong long-range interactions ^llj . The paradigmatic 
case of Newtonian gravity (i^ = 1) has served as a basis for developing methods for 
studying an important part of the phenomenology concerning the thermodynamics of 
systems with long-range interactions. In this regard, the isothermal spheres model 
has been widely used to study self-gravitating systems in the mean field limit [¥[|12|. 
Although in the mean field approach correlations are ignored, this model offers a 
mathematical tool for a suitable treatment of self-interactions in the system and turns 
out to be very accurate in the thermodynamic limit [T31II1], except near the critical 
points where the system collapses. It is also well known that self-gravitating systems 
possess equilibrium states with negative heat capacity, provided the system is isolated 
(microcanonical ensemble). Equilibrium in that case is ensured in a certain range 
in the space of parameters because isothermal spheres correspond to local maxima 
of the entropy with an extremely large lifetime that scales like the exponential of 
the number of particles: strictly speaking, they are metastable states [15]. In the 
microcanonical ensemble the system becomes unstable when the heat capacity passes 
from negative to positive, leading to the gravothermal catastrophe [2l[T6]. When one 
of these systems is put in contact with a heat bath (canonical ensemble), the range 
in the space of parameters where the heat capacity is negative is replaced by a phase 
transition [3J. In the canonical ensemble, mean field isothermal spheres correspond 
to states of local minima of the free energy and the isothermal collapse sets in when 
the heat capacity passes from positive to negative |12) . The self- gravitating gas has 
also been studied in the grand canonical ensemble with the mean field approach and 
Monte Carlo simulations. There, the instability sets in at a critical value of the 
parameter controlling the state of the system which is different from the critical values 
in the microcanonical and canonical ensembles |13l[T4] . This illustrates the fact that 
different ensemble representations are, in general, inequivalent and the thermodynamic 
behaviour of the system strongly depends on the control parameter used to specify its 
thermodynamic state. In addition, systems with an attractive interaction potential 
1/r'^ with < !/ < 3 in 3 dimensions were also considered in the microcanonical 
ensemble |17'fl8] and also there a critical energy was found bellow which these systems 
undergo a gravitational-like phase transition. 

Moreover, de Vega and Sanchez obtained [14] the equation of state of the self- 
gravitating gas and pointed out that it was customary to assume it without any 
derivation. To obtain the equation of state they took into account the condition of 
hydrostatic equilibrium and found that the system locally behaves as an ideal gas. The 
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same result is also found using the condition of hydrostatic equilibrium for a system 
with arbitrary long-range interactions in the mean field limit |19) . 

In this framework, our aim here is, on one hand, to analyze the local thermo- 
dynamics of systems with strong long-range interactions in d dimensions by computing 
the local thermodynamic potentials per particle in the mean field limit. By volume 
integration of the relation between the different thermodynamic potentials at the local 
level, we find the corresponding equation satisfied by the potentials at the global level. 
The local equation of state is obtained by computing the local entropy per particle 
and using the condition of local thermodynamic equilibrium instead of assuming 
hydrostatic equilibrium. Although in this case the two equilibrium conditions lead 
to the same result, both approaches are conceptually different. On the other hand, it 
is shown that the potential energy enters as a thermodynamic variable which modifies 
the global thermodynamic potentials. As a result, we find a generalized Gibbs-Duhem 
equation which relates the potential energy to the intensive variables of the system. 
For the marginal case where the power of the decaying interaction potential is equal 
to the dimension of the space, the usual Gibbs-Duhem equation is recovered. 

The rest of the paper is organized as follows. In section [21 we consider a 
system with strong long-range interactions in d dimensions and compute the local 
entropy, Helmholtz free energy and grand potential per particle in the microcanonical, 
canonical and grand canonical ensembles, respectively. We also consider the one- 
particle distribution function approach. In section[3J we obtain local relations satisfied 
by the thermodynamic potentials as well as the local equation of state. In section HI 
the global equation of state is obtained for a d-dimensional system and by integrating 
the local relations, an equation satisfied by the set of global thermodynamic potentials 
is found. We also obtain the generalized Gibbs-Duhem equation and consider some 
examples of its application. In connection with the latter, we point out a close 
relationship existing between our treatment and Hill's [20| thermodynamics of small 
systems. Finally, in section [Sj a discussion of our results is presented. We will use 
units such that fc^ = f . 

2. Mean field local thermodynamic potentials 

In this section we will compute the mean field local entropy, Helmholtz free energy 
and grand potential per particle in the microcanonical, canonical and grand canonical 
ensembles, respectively. Since the long-range character of the interactions depends on 
the dimension of the embedding space, we will consider the system is d-dimensional. 
Interactions are introduced through a long-range pair interaction potential (/'(<7i,<7j) 
which depends on qi and qj, the positions of particles i and j. Although here it is 
not necessary to make explicit the form of the interaction potential, we will formally 
assume that at large enough distances it behaves as 

Hqi,qj) ^ i^Wi- Qj]^" : (1) 

where k is a coupling constant and < v < d. This power law for the interaction, 
however, will be considered in Section 2] to derive global thermodynamic properties of 
the system. At the end of this section, we will briefly discuss how to compute local 
magnitudes with the one-particle distribution function approach, obtaining also the 
local entropy per particle. 
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2.1. Local entropy: microcanonical ensemble 

Consider a system of N classical point-like particles of equal mass m enclosed in a 
container of d-dimensional volume V ^ L'^, L being a characteristic length defining 
the size of the system. Positions and momenta of particles are respectively described 
by Qi = {q},qi, ■ ■ ■ , qf) and pi = {pl,Pi,. ■ ■ ,pf), 1 < i < N. The Hamiltonian of the 
system is 

HNiq^,p^)^Eo + W, (2) 
where the kinetic and potential energies, Eq and W, respectively, are given by 

^o=E^ and W^Y.<l>(l^,q,). (3) 

i i>j 

In the microcanonical description the state of the system is characterized by a fixed 
value of the energy E and the number of accessible microstates with energy E > Hn 
in d spatial dimensions is given by 

where d^'^^r = 0^=1 d'^Pid'^qi and 6{x) is the Heaviside step function. The domain 
of spatial integrations extend over the d-dimensional volume V and there are no 
restrictions on the domain of the momenta; this will always be the case throughout 
the paper unless another domain is specified. Accordingly, the microcanonical entropy 
is given by S{E) = \\\Y,{E). Introducing 

^ = ^^^' '^^^^ 
and integrating over momenta we have 

-\dN , 




N 



)dN/2 / 
^ \ £~e~'T.Hcl^,qJ) \ , (6) 

where d'^^q = YliLi '^'^Qi and r(a::) is the Gamma function. Using the integral 
representation [13] 

^ r(7+l) p exp(ifcx) 
x^O{x) = — / dfc . ^ , (7) 

^ becomes 



T.(E) ^ ^— / ^^-^ I dfc 
^ ' 2ttNI / ' 



— oo 



X exp 



iklS- 0(g,:, g,) ] - + 1 ) In(ifc) 



(8) 
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For convenience, let us introduce the parameter /3 according to e/3 — ifc and define 



1 P d'^^o / ^ dN \ 

^ J -^exp I -/3|]0(q„g,)-— ln(e/3) I , (9) 

so that ([5]) can be rewritten as 

^(E) r ^ exp [/3E + J(/5) - In . (10) 
J-ioo 27ri 

To evaluate ^ , we will employ a general method which was previously used in the 
context of long-range interacting systems to treat Newtonian gravity in d = 3 \^JT^ . 
The method assumes that for long-range potentials the interactions of particles with 
near neighbours are negligible as compared with the interactions with the rest of 
more distant particles. To implement this method, the d-dimensional volume V is 
divided in M cells, 1 < Af < N, of volume Va = V/M, so that J2aLi Va = V. 
Each cell is located at Xa = (a;^,^^, . . . ,xf) and contains Ua ^ 1 particles such that 
J2a=i ~ This construction is self-consistent if the number of particles is large 
enough, i.e., it is implicitly assumed that N oo. It is also assumed that the size of 
each cell is small enough so that the potential can be considered constant through the 
cell, but, of course, it varies from cell to cell. It follows that the total potential energy 
can be approximated as « J2a>b ''^a'nb<p{xa, xt) ■ 

Furthermore, with the help of the multinomial theorem it is not difficult to derive 
the following identity: 



{ni,...,nM i i=i *=1 

where {ni,...,nj\,f} means all possible values ni, n2, . . . , riM of the occupation 
numbers and the Kronecker delta restricts the total number of particles to TV. 
Equation pT|) can be seen as a way to perform an integral in the dTV-dimensional 
configuration space by summing over all possible occupation number distributions 
through the cells of a single-particle configuration space. Thus, using the multinomial 
expansion one approximates ([9)) by 



eJ(/3) 



{ni,...,nM} i=i «— J- 



X exp -/3 y nanb(t>{xa,Xb) - - ln(£/3)Vn, . (12) 





The above approximation formulated in a discrete way allows one to obtain a 
continuous field representation. This continuum limit is obtained by introducing the 
number density, whose value in each cell is given by n{x) = Ua/va- Consequently, 
summations over the cells become spatial integrals and the summation over occupation 
number distributions becomes a functional integration in the number density: 

— and J2 ^ ■ (^3) 

a=l ° {ni,...,nM} 
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Using these prescriptions the total potential energy in the mean field limit takes the 
form 

W[n\ = ^ y" n{x)n{x')<j){x, x') d'^x dV . (14) 
In addition, with Stirling's approximation we also have 

M 



expl- / n(a;) In d'^a; ) . (15) 



TT-f- 

-IJ- n,! \V 

i—l 

Therefore, taking into account that the Kronecker delta restricting the number 
of particles becomes a Dirac delta in the continuum limit and using the integral 
representation 5 {x) — (27ri)^^ P^oo dae~"^, one gets 



e^(/3) 



lOO 

ioo 

da 



27ri / 2^"- exp (^J[?i;a,/3]) , (16) 
where we have introduced the functional J[n; a, /3] defined by 
J|„;o,^l.a(A'-/..WA) -/„,.),„ !l<HlId'. 

j n{x)n{x')(l){x,x')d'^xd'^x' -'^\n{el3) j n{x)d'^x. (17) 
Hence, by replacing ([TBI) ([TOll - we have 

Vn exp {S[n]a,ri\^ , (18) 



/"IOO 

^S(E) _ 



da r°° d^ 
27ri ./ 27ri 



where we now have introduced the dimensionless variable rj = eji and 

S[n\ a, rj] = rje'^E + J[n; a, yye"^] =/?£' + J[n; a, P] . (19) 

In the definition of S we have neglected the term In (e/3) = ln(r/) coming from the 
exponential in (jlOl) . since it does not contribute in the limit — > oo. This is the same 
as assuming that the microcanonical entropy is given by the logarithm of the number 
of microstates or. equivalently, by the logarithm of the density of states. To proceed 
further, p8p will be evaluated using the saddle point approximation. That is, the value 
of the integral will be given by the exponential of S[ns; as,ris], where ns{x), as and 
77s are the number density and the value of the parameters that maximize S[n;a,r]]. 
Accordingly, the microcanonical entropy can be approximated by S{E) w S[ns; as, r?s]- 
In order to simplify the notation, in what follows we will omit the subscript s, since we 
will only consider the stationary solutions given by the saddle point approximation. 

The parameters a and 77 in (llSp can be viewed as Lagrange multipliers that 
respectively restrict the value of the number of particles and the energy. Indeed, 
the mean field equations in the saddle point approximation given by the extremal 
condition ((5„; da, drj) S[n; a, 77] = (0; 0, 0) yield 



N = J n{x) d'^x, (20) 
J n{x) d'^x + W[n] (21) 



2(3 



after computing the derivatives of S[n; a, rj] with respect to a and rj. Moreover, taking 
into account (jl7p . variations of the number density in (|19l) lead to 6S[n\ = 6J[n\ with 

(5 J[7i] = J -I3^{x) - a - In (^72{x)Vt]'^/'^^ Sn{x) d'^x , (22) 
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where the term containing the potential ^(x) comes from the variation 5M^[n] = 
J ^{x)Sn{x) d'^x and is given by 

^{x) [ n{x')(f>{x,x')d''x' . (23) 



In terms of $(a;), the total potential energy takes the form 

W[n] = ^ y" n{x)<i>{x) d'^x . (24) 

Hence, by setting the functional derivative with respect to n{x) to zero, (5„S'[n] = 
J„J[n] = 0, it follows the coupling between the number density and the interaction 
potential through the relation 

n{x) = Aj;'' exp [-^ {<P{x) - //)] , (25) 

where we have introduced At — V^l'^r\^l'^ and /i — —a/ (3. As a result, by replacing 
(|20p and (PT|) in ([T7]) and (|19l) , the extrema of S and consequently the microcanonical 
mean field entropy become 

2 + d^ 



S = j n{x) -ln(n(a;)A^) + 



d^a; , (26) 



which is an integral over the volume and therefore the integrand can be interpreted 
as the density of entropy. From this density of entropy one can construct the local 
entropy per particle and in this way, as we shall see below, information about the local 
thermodynamic nature of the system can be obtained. Taking into account (j26p . the 
local entropy per particle is given by 

s(a;)--ln(n(a;)A^) + ^, (27) 

in such a way that the total mean field entropy can be written as 

S = j n{x)s{x)d'^x. (28) 

In addition, using p7)) and p^ . it can be checked that /3 is the inverse 
temperature, /i is the chemical potential, 

1 fdS\ „ n fdS\ 

and that 



is the thermal wavelength. 

The microcanonical mean field entropy ([26| is an extremal solution of the mean 
field equations. In order to guarantee thermodynamic equilibrium, we will assume 
that there exist a certain range of parameters where S" is a maximum and thus, our 
analysis here must be understood restricted to such a range. It is not difficult to see 
that such a range in the space of parameters does exist: if the temperature is high 
enough so that /3$(a;) <C 1, the number density becomes n w N/V and the stable 
global ideal gas behaviour is recovered. The critical values of the parameters that set 
the instability will depend on the interactions under consideration. 
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2.2. Local Helmholtz free energy: canonical ensemble 

Here we will consider that the control parameter used to specify the equilibrium 
configurations is the inverse temperature /3. At fixed temperature one must use the 
canonical description and in this section we will compute the canonical partition 
function, Z{/3), in the saddle point approximation. In doing so, the canonical 
Helmholtz free energy, A(/3), can be obtained and written as an integral over the 
volume which allows us to identify the local free energy. 

The canonical partition function can be obtained from the microcanonical density 
of states. We = ds^iE), by computing its Laplace transform. The density of states 
is then given by the inverse Laplace transform of the partition function, 



.S(i?)= r ^e'^i^e^^. (32) 
J-ioc 27ri 



By computing the derivative of expression ()10[) with respect to E one gets 
dE' 

Therefore, the comparison of both expressions (|3T]) and p2l) leads to the identification 

Z{I3) = e-"^^^ , (33) 

and thus, one sees that the function J(/3) defined in © is just the Massieu function 
related to the Helmholtz free energy via J(/3) = — /3A(/3). Hence, the major part 
of the work has already been done writing down the microcanonical description, 
since the canonical partition function is given by It is clear that the extremal 

condition ((5„; 9q) J[n; a, /3] = (0;0) allows one to get an approximate expression 
for the Massieu function as J(/3) ~ j[n;a,/3], where, as before, n{x) and a are 
given by the saddle point equations. The saddle point equation associated with 
the parameter a constrains the number of particles leading to the condition ([^0]) . 
Of course, in the canonical description only the temperature is fixed and there is 
no second Lagrange multiplier enforcing the condition ([?!]). Instead, according to 
the saddle point approximation and using (jl7p . one sees that the mean value of the 
Hamiltonian {E) = —dp In Z w —dp J = E, becomes 

E=— n{x) d'^x + W[n] (34) 



2/3 

in the mean field limit. 

On the other hand, in order to obtain the relation between the number density 
and the interaction potential, variations with respect to the number density of J[n\ 
has to be computed. This task has also been done in the previous section and 6J[n] is 
given by expression As a result, the distribution of particles that maximizes the 
microcanonical entropy also maximizes the Massieu function in the canonical ensemble 
and the number density is then given by Hence, using P7)) . (|^D|) and the 

saddle point approximation allows one to write the canonical mean field Helmholtz 
free energy A — —TJ as an integral over the volume in the mean field limit: 



X . (35) 



A^ J n{x) |r [In {n{x)X^) - l] + ^$(a;)| d'' 
Therefore, the mean field local Helmholtz free energy is readily identified, 

a{x)^T[\n{n{x)X^) -l]+^^{x), (36) 
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since 

j n{x)a{x)d'^x. (37) 

The canonical mean field entropy can be obtained via a Legendre transformation 
of the canonical mean field free energy, S — l3{EA), and, as a consequence of the 
saddle point approximation [5] , it coincides with the microcanonical mean field entropy 
(|26p. The mean field solution in the canonical ensemble is valid only in the range of 
parameters where J is maximum. 



2.3. Local grand potential: grand canonical ensemble 

Consider the canonical partition function as a function of the number of particles and 
its Laplace transform: 



/>oo 

C[Z]{a) = / dNZ{N)i 
Jo 



(38) 



Taking TV as a continuous variable, the canonical partition function is constant with 
value Z{N) for all N in the interval [N, N + 1). Therefore, the integral in ((38|) can be 
written as a sum of integrals over the intervals [N, iV + 1] and subsequently integrated: 

C[Z]{a) ^ ^(^) / dTVe""^ (39) 

N=0 

1-e- 



N 



a 

N=0 



J2 e~"^Z(7V) . (40) 



From here, since a = —PfJ-, one recognizes the grand canonical partition function 

oo 

Z(a) = 5]c-"^Z(7V), (41) 

which we will take as the starting point to derive the thermodynamics of the system. 
Introducing C(a) = [1 — exp(— a)]/a and applying the inverse Laplace transform to 
(gni) one gets 

/ioo 1 
^C(a)Z(a)e"^. (42) 
-ioo 27ri 

Since in the previous section we obtained 

/ioo A f / \ 

— J Vnexp(^J[n;a,(3]j , (43) 

the grand canonical partition function reads 

Z = Jvn exp (j[n; a, /3] - aiV - InC(a)) . (44) 

The term In (^(a) in the exponential above, which is a correction due to the discreteness 
of A'', can be safely neglected in the limit N ^ oo. As before, we shall evaluate the 
integral (|44|) using the saddle point approximation, hence Z « exp( J — aN). In this 
case, only the number density such that Sn{J[n\ — aN) = SnJ[n\ = is required. Thus, 
one more time the number density is given by Furthermore, the mean value of 
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the number of particles can be written as {N) — —da InZ w —da{j — aN) = N and 
using ([T7)) one gets 



N = 



J n{x) d'^x . (45) 

Analogously, the mean value of the energy takes the form {E) — — 9^1nZ w 
—dp{J ~ aN) ~ E, where, as in the canonical case, E is given by Introducing 
the grand potential 

n^-T\nZ, (46) 

the saddle point approximation yields — /3J7 « J — aN, which by using (fT7)) and ([M)) 
leads to 

n = J n{x) |r [In {n{x)X^) - l] - /.t + ^*(a;)| d'^x (47) 

in the mean field limit. Consequently, the mean field local grand potential per particle 
reads 

uj{x) ^ T [In (n(a;)Af,) - l] - m + \^{x) , (48) 



n ^ n{x)uj{x)d'^x. (49) 



such that 



The grand canonical mean field Helmholtz free energy is obtained by a Legendre 
transformation of the form A = H + fiN, which turns out to be equal to the mean field 
Helmholtz free energy in the canonical ensemble (|35l) . With a subsequent Legendre 
transformation, we would obtain the same mean field entropy as in the microcanonical 
ensemble S = /3{E — ^N — fi). As mentioned before, this is a consequence of the 
saddle point approximation. Although each ensemble has its own range of validity in 
the corresponding space of parameters, the thermodynamic potentials computed in 
the different ensembles have all the same form in the mean field limit. In the grand 
canonical ensemble the mean field solution is valid only in the range of parameters 
such that J — aN is maximum. 

2-4-. One-particle distribution Junction approach 

Let us consider a system of N classical particles of mass m confined in a container of 
volume V. These particles interact via a pair potential 4){x,x') which describes the 
interaction between particles at positions x and x' . The distribution of particles is 
given by the stationary distribution function / = f{x,p) defined in the 2(i-dimensional 
single-particle phase space, p being the momentum. The spatial distribution of 
particles is given by 

n{x) = J f{x,p)d''p (50) 

and the distribution function is normalized such that 

N = J f{x,p)d'^xd''p = J n{x)d'^x. (51) 

The energy of the system can be written in terms of the distribution function according 
to 

E = Eo + W = J + ia>(a;)^ f{x,p) d^x d^p , (52) 
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where the potential reads 

^{x) = J (j}{x,x') /(a;',p')dVdV (53) 
and therefore, the total potential energy takes the form 

W =^ J $(a;) f{x, p) A'^x A'^p (54) 

= iy c^{x,x')f{x,p) /(a;',p')d''a;d'^pdVdV. (55) 

In view of ([501) . it is clear that the potential and the total potential energy are given 
by equations and ([M)) . respectively. 

Assuming that the thermodynamics of the system is obtained from the Boltzmann 
entropy defined by 



S^- j f{x,p)\nl^^d''xd''p, 



(56) 



the stable configurations are those having a distribution function that maximizes this 
entropy. Here /c is a constant which does not depend on the thermodynamic state 
of the system. The distribution function which makes the entropy be extremal is 
obtained by solving the variational problem 

6S[f] ~ a5N[f] - l36E[f] = , (57) 

where a and /? are Lagrange multipliers enforcing the constraints ([FT]) and ([5^ . 
respectively. The Lagrange multipliers, as before, can be identified with the chemical 
potential, /i — —a//3, and the inverse temperature, /3 — T~^. One finds that the 
solution to the variational problem is the Maxwell-Boltzmann distribution given by 
/ = /cexp {-[1 + a + /3$ + /3|p|V(2m)]}, which can be rewritten as 



f{x,p)^Xrj.''[- exp 



13 



d/2 



^/?|^ -/?(<!>(.)-,) 



(58) 



where = i^/ fcY^'^ /V^irmT. Integration of the distribution function ([58|) over 
momenta gives the local distribution of particles in the spatial domain, 

n{x)^ y/(a;,p)d'^p = A^''exp[-/3(<i>(a;)-/.)] , (59) 

which coincides with ([251) if oi^^ set the constant /c in Xt to be — e/{2TihY. One 
also sees that 

"'^ r lr,|21 

(60) 



f{x,p) = n{x) ( ) exp 



27rm 



2m 



As in the case of the density above, any local magnitude can be obtained by 
integration in phase space of the corresponding microscopic contribution weighted 
down by the distribution function. That is, any local magnitude per particle q{x) 
associated with the global quantity Q can be defined through 



provided 



nix)q{x) = J fix,p)Qix,p)d''p (61) 
Q = J f{x,p)Q{x,p)d''xd^p. (62) 
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In the case of the entropy, e.g., one takes S{x,p) — — ln(/(a;,p)//c) which after 
integration using (pD|) gives 



J f{x,p)S{x,p) d'^p^n{x) 



In {n{x)X'^) + ^ 



(63) 

The term in brackets coincides with (|27|) . the mean field local entropy obtained in the 
ensemble representations. By Legendre transforming the mean field entropy, all the 
other mean field thermodynamic potentials can be obtained. Thus, the ensemble 
representations and the one-particle distribution function approach give the same 
mean field local thermodynamic potentials. 

Depending on the boundary conditions given by the control parameter, the 
thermodynamic equilibrium is ensured if the corresponding thermodynamic potential 
is maximum. If the system is isolated with fixed energy E, then the conditions 
on the stability of the mean field solution are found by considering second order 
perturbations of the entropy S. On the other hand, if the system is in contact with 
a heat bath which fixes the temperature, second order perturbations of the Massieu 
function J = S—(3E must be considered, with E = E given by ([52| . If, in addition, the 
number of particles fluctuates constrained by a given value of the chemical potential, 
the stability is guaranteed if the quantity J — aiV, with N — N given by ([5T|) . is a 
maximum. Therefore, in this case second order perturbations of J — aN have to be 
taken into account. 



3. Local relations and the local equation of state 



In the previous sections we wrote the thermodynamic potentials as integrals over 
the volume of the system, which leads to a natural definition of local quantities per 
particle. Analogously, using (|21l) and (|24|) . the local energy per particle takes the form 



e(a;) = eo + i$(a;) , 
where eo is the local kinetic energy per particle. 



eo 



d 

-T . 

2 ' 



such that the total energy and kinetic energy are given by 

n{x)eo d'^x . 



We also introduce the local volume per particle, v{x), given by 

1 



v{x) = 



(64) 

(65) 

(66) 
(67) 

(68) 



n{x) ' 

since, obviously, V = J n{x)v{x) d'^x. 

Once local variables are defined, the local equation of state can be obtained. In 
order to do that, it is useful to write the local entropy per particle in a more convenient 
way. Using ^U^, and (UHl), from (gT)) one obtains 



s{x) = In 



d/2' 



2 + d 



(69) 
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which is the Sackur- Tetrode equation in d dimensions formulated in terms of local 
variables. The local entropy per particle is therefore explicitly obtained as a function 
of the local variables Cq and v, s = s(eo, v). Thus, in the equilibrium framework, one 
infers that the local internal energy only has contributions coming from kinetic degrees 
of freedom while self-interactions play the role of an external field that perturbs the 
gas. In addition, local thermodynamic equilibrium implies the relations |2H[22] 



where the second of these relations is the local equation of state of the system. In our 
case, this local equation of state is indeed the one corresponding to an ideal gas, 

p{x) = n{x)T, (71) 

which is valid for any long-range pair interaction potential that can be suitably 
represented with the mean field approach. Since the mean field local entropy per 
particle takes the same functional form in the microcanonical, canonical and grand 
canonical ensembles, the above equation of state is valid in the three ensembles. Notice 
that the condition of local thermodynamic equilibrium together with expressions (|70| 
are formulated as a hypothesis in the framework of non-equilibrium thermodynamics, 
but here they are trivially satisfied provided the whole system is in equilibrium. 

To get more insight in the relation between local thermodynamic variables, we 
use ((25)) and write the chemical potential as 

^i = lJLa{x) ^ <^{x) , (72) 

where 

^if){x) =T\n{\'^n{x)) (73) 

possesses the same functional dependence on the number density as the chemical 
potential of and ideal gas, but with the density given by (|25|) . Also using ((25| . the 
local entropy per particle can be rewritten as 

s{x) = ^{^T + ^{x)^+l-^^, (74) 

and because p{x)v{x) — T one obtains 

Ts{x)^[eo + <^{x)]+p{x)v{x)- ^i. (75) 

The term in brackets is the total energy of a particle at the point x\ locally the 
potential ^{x) acts like an external field. However, the difference between $(a;) and 
an authentic external field becomes manifest when one sums the contribution of the 
whole system: multiplying [eo -I- *I'(a;)] by the density and integrating over the volume 
does not give the total energy. In order to obtain the total energy one must take 
into account that the total potential energy is due to self-interactions and this is 
the reason why one defines the local energy per particle according to ([64]) (the total 
potential energy is a functional quadratic in the density). In this respect, using (|64p . 
equation (|75p can be alternatively written in the form 

Ts{x) ^ e{x) + p{x)v{x) - ^ + i$(a;) , (76) 

where the last term highlights the fact that self-interactions are actually considered. 
Nonetheless, the statement that particles behave locally as an ideal gas can also be 
seen by rewriting (|75p explicitly in the form 

Ts{x) ^ eo+p{x)v{x) - iio{x) . {11) 
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The local Helmholtz free energy per particle, a{x), can be obtained by the 
Legendre transformation 

a{x) = e{x) - Ts{x) = ^ - p{x)v{x) - ]^^{x) . (78) 



It can be checked that (|78)) coincides with ((36)) . The local Helmholtz free energy per 
particle can also be written in the form 

a(a;) = ao(a;) + i$(a;) , (79) 

where ao{x) is the local Helmholtz free energy per particle of an ideal gas at the point 
x, which is given by 

ao(a;) = -T [in {v{x)X;j.'^) + l] . (80) 

Hence, it also satisfies 

ao{x) = eo - Ts{x) = ^q{x) - p{x)v{x) . (81) 

Moreover, the local grand potential ((48|) is equivalent to the Legendre transformation 
ui{x) = a{x) — and with elementary manipulations it can be rewritten as 

Lu{x) ^ ~p{x)v{x) - ]^'^{x) . (82) 

The local grand potential per particle of an ideal gas takes the form 

u}q{x) = ao{x) - noix) = ~p{x)v{x) , (83) 
so that it is related with u!{x) through 

cj{x) ^ ujo{x) - ^^x) . (84) 

By replacing (|82p in and using the local equation of state one sees that the total 
grand potential is closely related to the kinetic and potential energies, 

n = - J p{x)d'^x-W ^-^Eo-W . (85) 

It is worth noting that the Legendre transformation of a local variable with 
respect to the local volume, in general, does not agree with the corresponding Legendre 
transformation of the global variable with respect to the total volume. The reason 
is that, in general, the pressure is not uniform, in contrast with the other intensive 
variables /i and T. For instance, integration of the quantity 

g{x) EE a{x) + p{x)v{x) ^ fi - ^<^>{x) , (86) 

that is, 

G = y n{x)g{x) d'^x = A+ "^Eq = fiN - W , (87) 

does not coincide with the total Gibbs free energy G = A + PV when long-range 
interactions make the system be inhomogeneous. Here P — p{x)\x ^ov is the pressure 
evaluated at the boundary of the d-dimensional system. The quantity PV depends on 
the total value of the potential energy and its functional relation with the total kinetic 
energy; e.g., in the case of Newtonian gravity in d = 3 one has 2Eo + W — 3PV. An 
analogous situation occurs if the enthalpy H = E + PV is considered. 

In the next section, the global equation of state and an expression for the total 
Gibbs free energy will be obtained for systems in spherical containers with interaction 
potentials of the form ([T]). In such cases, the local Gibbs free energy can be trivially 
derived from the global free energy G. 
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4. Global magnitudes and the generalized Gibbs-Duhem equation 

Useful relations between global quantities can be obtained by integration of the local 
equations. Multiplying both sides of (|76[) by n{x) and integrating over the volume 
yields 

TS = E+-,Eo-tiN + W. (88) 
d 

However, it would be desirable that the quantity PV appears explicitly related to the 
other thermodynamic quantities. In order to do that, the form of the pair interaction 
potential has to be considered. We shall see that once the quantity PV is explicitly 
introduced, not only the Gibbs free energy can be computed but also the Gibbs-Duhem 
equation satisfied by the system. 

We will consider that the container is a (d — l)-sphere of radius R, in such a way 
that particles are confined by the boundary condition (^qj) + (g,^) + • • • + [qf) < i?^, 
1 < i < N. The 0-sphere is the pair of endpoints of the line segment of length 2R. 
Besides, the volume of the system reads V = 7r'^/2i?'*[r(d/2 + 1)]"^ 

Thus, from using coordinates defined according to — q-jR and rescaling 
energies by introducing P^[T5] 

^ ^ _ _^ / r(d/2 + i) \ 

i?" In 

^^"■^'""^■^ = pv^'^^^^' ^■'■^ " iv^ H ~ ' ^^^^ 

the entropy takes the form 

YN(2-u)/2 

S^\n +lnE(A). (91) 

Here 

"^"U'l^l^y lr(d/2 + i)j 

and 

\'~/' / A. \ 

A-^<^(r„r,) WA-^<^(r„r,) . (93) 

i>3 ) \ i>3 / 

Notice that E(A) depends on the volume only through A, since in the ri-coordinates, 
which belong to a certain domain D, the limits of the integrals in (|93p are pure 
numbers. Taking into account that {dE^)v = ^/E and (SyA)^; = i'A/{dV), one has 

1 fdS\ A91nS(A) 

- ' ' - (94) 



T \dE},, E dA 
P 
T 



fdS\ N iy\ z.A91nS(A) 

where P — p{x)\x ggy. From the above equations and using that E = Eq + W one 
gets 

PV W 

NT = '^'dNf^ ^''^ 
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which is the exact microcanonical equation of state of the system. However, in the 
mean field limit this equation holds also in the canonical and grand canonical ensembles 
if mean values are taken such that E = E and N = N. Of course, the last statement is 
only valid in the domain of the space of parameters where each ensemble representation 
is well defined. With this in mind, for simplicity, in what follows we will not distinguish 
between E and E or N and N. 

From the equation of state (|96|) the virial theorem is obtained. 



2Eq + vW ^ dPV , (97) 

which is particularly useful to express the relation between global quantities. It is also 
useful to introduce the long-range parameter a defined as 

d — v , , 

a = — — , < cr < 1 , 98) 
d 

which together with (|97p enables us to rewrite (|88p in the form 

TS = E + PV - fiN + aW . (99) 

The marginal case v = d correspond to systems with long-range parameter cr = 0, 
so that the above equation and the thermodynamic potentials obtained below reduce 
to the usual ones in short-range interactions thermodynamics. The Helmholtz free 
energy, the grand potential and the Gibbs free energy are readily obtained from 
and take the form 

A = E - TS = ^^N - PV ~ aW , (100) 
n = A - iiN ^ -PV - aW , (101) 
G = A + PV ^ ^iN ^ aW . (102) 

The above expressions generalize previous results for the self-gravitating gas [T51ll4| . 
In addition, by differentiating (|M)) one gets 

TdS = dE + PdV - ^^lAN + adW - iVd^ - SdT + VdP (103) 

and since TdS = dE + PdV — iidN we must have 

adW = 5dr - VdP + Ndfi , (104) 

which is the generalized Gibbs-Duhem equation satisfied by the long-range interacting 
systems we are considering here. For systems with long-range parameter a — 0, i.e. 
the marginal case v = d^ the usual Gibbs-Duhem equation is recovered. From (|104p 
we infer that when long-range interactions are present, the intensive variables are 
independent. Moreover, in view of (|104p . the following relations are obtained for 

'(f).-; "(fL-- "(f), ^ 

The above expressions imply, for instance, that the entropy can be obtained from a 
derivative of the potential energy. Since such thermodynamic relations have not been 
reported previously for systems with long-range interactions, it would be interesting 
to check them by considering some solvable examples. Next, we will consider the case 
where the interaction potential is spatially uniform and the case where the system 
is a self- gravitating gas. After these two examples we will expose a common feature 
between systems with long-range interactions and small systems. 



N . (105) 
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Spatially uniform interaction potentials in d dimensions 

These systems are characterized hy v — such that the system is homogeneous and 
(7 = 1. We simply have cj) — k, hence $ — kN and W — kN'^/2. In order to apply 
(jlOSp . the potential energy must be written as a function of the intensive variables, 
W = W{T, P, n). Using ([25| and the equation of state (|96l) . which in this case reads 
PV = NT, one obtains n = N/V = P/T = X^'^ exp [-^{kN - fi)] and the number of 
particles becomes 



N : 



Tin 



The potential energy then takes the form 

wiT,p,^l) 



1 

2k. 



Tin 



(106) 



(107) 



Taking into account that the local entropy per particle f27|) is uniform in this case, 
whence S — Ns, it is straightforward to see that 



S 



Tin 



The remaining two relations p05p follow from V = NT/P and (|106p . 



4-2. Self- gravitating isothermal spheres in d — S 

The calculation of global thermodynamic quantities of self-gravitating isothermal 
spheres is well understood and here we will just consider it as an example to illustrate 
the relations (|105p . In particular, the first of these relations will be used to obtain the 
entropy of the self-gravitating gas. We will write down only the necessary expressions 
to carry out our task and refer the reader to [4l[T2j[T5] for details. 

In the self-gravitating gas we have — 1, thus a = 2/3, and k = —GN'm'^, where 
Gat is Newton's constant. Also, Newtonian systems satisfy the Poisson-Boltzmann 



equation which, 

and ip = f3 {^{x) 



after introducing the dimensionless variables ^ — (47r|K|/3 n(0)) 
— $(0)), becomes the Emden equation 



1/2, 



(108) 



Here x — \x\ and /3 = 1/T while n(0) and $(0) are the density and the potential 
at the origin. Global thermodynamic magnitudes are expressed in terms of these 
functions evaluated at the boundary of the spherical container, thus, for convenience 
one introduces = (47r|K|/3 n(0))^/^i?, -00 = V'(^o) and ip^ = '0'(?o): where primes 
denote the derivative with respect to ^. Moreover, with a suitable change of variables 
the Emden equation can be transformed into a first order differential equation [23j. 
When such variables, usually denoted as (w, u), are evaluated at ^ = and "0 = V'Oj 
with ip' = 001 they read 



VO = ^000 



and 



Uq 



and therefore satisfy 

duo 
dvo 



vo{uo - 1) 



(109) 



(110) 
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On the other hand, it can be shown [Hll2| that the inverse temperature and the 
potential energy can be written as 

1 R , , 

(111) 



T \k\N 

W = NT{uo-3) . (112) 

What we need is to express the potential energy as a function of T, P and /i only. 
Since $(i?) = kN/R, equations ^ and ^ yield 



R 



Tin 



PAf, 
T 



(113) 



where we have taken into account the local equation of state so that P = n{R)T. 
From pi3p and (lllip . one sees that 



vo(T,P,fj,) = In 



and therefore 



I dT 



/i 



T 



T \ T 



(114) 



(115) 



Combining ()112p and the global equation of state (j96| one obtains 3PV = NTuq and 
hence, taking it into account and using (jllip to express N , the radius can be written 
as 

/ \ 1/2 

^ ^ . (116) 



(117) 



i? = T , 

V47r|K|P, 

Using (|llip and (|116p . the potential energy (|112p takes the form 



WiT,P,^) 



( 3/2 



3ui 



1/2 



(47r|K|3p)^ 

which depends only on the desired variables, since uq = {vo{T, P, /i)). We then have 



\dT J 



T 2va \dT 



1 



vojun - 1) dup 
uo{uq - 3) dvQ 



(118) 



Therefore, according to ()105p and using ()110p . (|112p and (|115p . from the above 
equation (jllSp the mean field entropy is obtained. 



S = -N 



- - 2uo + - 
T " 2 



(119) 



1 /2 

Using that the thermal wavelength can be expressed as At — \2iTh-^ Rvo/ {N\K\m)\ 
and that P/T = Nu(i/{W), from ITTi)) one gets 

?3™,3 I 1 3 ^ 



1 /27V^W|^N 

2 1, j ' 



and thus the entropy becomes 



vq + 2uo - - In(wo) - In(woWo) - 3 



•S'o ) 



(120) 



(121) 



as given in [Ti|, with 5o = iV/21n [2NR^m^\K\^l (Trefi^)]. This verifies the first of 
the relations (|105p for the self-gravitating gas, and the remaining two follow from an 
analogous procedure. 



Local thermodynamics and the generalized Gibbs-Duhem equation 19 
4.3. Small systems 

Systems with a small number of constituents can be treated from a thermodynamic 
point of view by considering them as independent members of an ensemble. Since such 
an ensemble of independent small systems is a macroscopic system itself where the 
standard thermodynamics applies, thermodynamic properties of a single small system 
can be derived. This is the well known approach introduced by Hill [20]. The basic 
idea is to consider a 'chemical potential', denoted as S", that accounts for the energy 
gained by the system when the number of members of the ensemble, ^ , varies. The 
energy, volume, number of particles and entropy of each small system are E, V, N and 
S, respectively, while the total magnitudes corresponding to the whole ensemble are 
given by Et = ^E, Vt = J^V , Nt = JVN and St = ^S. Since the whole ensemble 



is a macroscopic system, these magnitudes satisfy 

TdSt = dEt + PdVt - ^idNt - Sd.yV (122) 

and hence, TSt ~ Et + PVt — fiNt — ,yK. Thus, for any small system: 

TS = E + PV - ^iN - S" , (123) 

TdS = dE + PdV - ^idN , (124) 

d,^ ^ -SdT +VdP - Nd^. (125) 



Hill mentioned that the use of different environmental variables, i.e. control 
parameters, would lead to different descriptions of the thermodynamic phenomena 
when small systems are considered [20] . In this case, ensemble inequivalence is due 
to the finite size of the systems. In fact, negative heat capacity leading to ensemble 
inequivalence is a generic feature of finite systems at phase coexistence [6J. On the 
other hand, systems with long-range interactions can be regarded as finite systems as 
well, in the sense that the range of the interaction is comparable with the size of the 
system [B]. 

Our aim here is to point out another common feature between small systems 
and systems with long-range interactions. Due to non-negligible self-interactions, an 
ensemble of independent systems with long-range interactions cannot be considered. 
Nonetheless, if one sets ^ = 1, fixed, so that d^y — 0, and relates the energy S' with 
the potential energy such that (0 = —aW, equations (|123p and (|125p become the same 
equations as ((99)) and (|104p . Although W and S" are introduced for different reasons, 
they formally play the same role. Therefore, a small system such as a macromolecule 
and a system with long-range interactions such as a star, both being finite system, 
can be formally described by the same thermodynamic relations. 

5. Discussion 

We have studied the local thermodynamics of systems with strong long-range 
interactions in d dimensions by computing local thermodynamic potentials per particle 
in the mean field limit. We have considered pair interaction potentials that decay as 
l/r", with < v < d. By computing the local entropy per particle and using the 
condition of local thermodynamic equilibrium we have obtained the local equation of 
state, which corresponds to the isothermal ideal gas equation, regardless of the long- 
range interaction potential that couples the constituents of the system. This result 
coincides with the results that were obtained in previous works where the condition of 
hydrostatic equilibrium was assumed [141119] . Thus, we have shown that interactions 
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locally play the role of an external field that perturbs the gas and that the local 
entropy per particle is a function of the local kinetic energy and local volume: it is 
the Sackur- Tetrode entropy formulated in terms of local variables. The local entropy 
per particle was computed in the microcanonical ensemble and with the one-particle 
distribution function approach, leading to the same result. 

Different ensemble representations possess, in general, different ranges of validity 
in the space of parameters specifying the state of the system. Since any global 
thermodynamic potential computed in different ensemble representations has the same 
functional form in the mean field limit, local thermodynamic potentials have also the 
same functional form in the different ensemble representations. By volume integration 
of the relation satisfied by the thermodynamic magnitudes at the local level, we 
have obtained the equation satisfied by the set of global magnitudes. Remarkably, 
the potential energy enters as a thermodynamic variable which modifies the global 
thermodynamic potentials. In the thermodynamics of systems with short-range 
interactions the potential energy is obviously involved since it is included in the 
total energy, but it does not appear as an additional contribution which modifies 
the thermodynamic potentials. Such a contribution is proportional to a = (d — i/) jv 
and hence, for the marginal case p = d, the usual thermodynamic relations valid for 
short-range potentials are recovered. 

As a consequence of this contribution coming from the total potential energy, the 
Gibbs-Duhem equation is modified and the intensive variables can no longer be seen as 
independent if long-range interactions are present in the system. By taking advantage 
of the thermodynamic relations obtained from this generalized Gibbs-Duhem equation, 
the entropies of a system with a spatially uniform potential and of the self-gravitating 
gas have been computed. 
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